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ABSTRACT 

The  vortex  valve  pneumatic  oscillator  is  described  .  Data 
obtained  for  variations  of  inlet  area,  bias  flow  and  temperature 
are  presented.  The  input  impedance  of  the  valve  is  determined 
experimentally  and  found  to  have  a  region  of  negative  resistance. 
If  the  valve  is  supplied  through  a  capacitance  (a  volume)  the 
combination  produces  a  limit  cycle  oscillation  where  pressure  in 
the  volume  is  a  function  of  time. 

Theory  predicts  the  frequency  of  oscillation  to  be  insensi- 
tive to  bhe  supply  pressure.  Analysis  of  the  supply  pressure 
versus  frequency  data  shows  the  oscillator  is  nearly  pressure 
insensitive • 

Theory  predicts  the  frequency  of  oscillation  to  be  directly 
proportional  to  temperature  *     The  static  valve  characteristics 
were  found  to  vary  with  temperature.  Therefore,  the  temperature 
sensitivity  of  the  oscillator  was  predicted  from  the  static 
valve  characteristics.  The  ratio  of  predicted  temperature  sensi- 
tivity to  actual  temperature  sensitivity  was  3*2  to  1. 
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CHAPTER  I 
INTRODUCTION  AND  OBJECTIVES 

1.1  Introduction 

The  vortex  valve  is  a  pure  fluid  modulator  or  amplifier.  Tho 
operation  of  the  device  depends  on  the  dynamic  interactions  betv;een 
streams  of  fluid  rather  than  moving  parts.  The  absence  of  moving 
parts  and  the  absence  of  interface  problems  make  the  pneumatic 
oscillator  attractive • 

However,  the  sacrifice  made  for  having  no  moving  mechanical 
parts  is  difficulty  in  designing  the  geometric  shapes  analytic- 
ally. In  fact,  an  analytical  procedure  for  complete  design  does 
not  exist.  Therefore,  a  great  deal  of  the  design  is  accomplished 
emperically. 

1.2  Objective 

The  objective  of  this  thesis  is  to  design,  construct  and  test 
a  pneumatic  oscillator  whose  frequency  of  oscillation  is  insensi- 
tive to  supply  pressure  but  sensitive  to  air  temperature .  Design- 
ing a  device  which  is  pressure  insensitive  eliminates  the  need  for 
regulation  of  the  supply  pressure.  A  fluidic  temperature  sensor 
is  highly  desirable  for  temperature  sensing  in  systems  operating 
with  oxidizing  gases  or  gases  at  extremely  high  temperatures. 
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An  ideal,  application  for  such  a  device  would  be  as  a  c     ent 
of  a  fluidic  speed  control  system  for  a  gas  turbine. 

The  plan  of  attack  was  straightforward!  design  an  oscillator 
that  is  theoretically  pressure  insensitive  and  temperature  sensitive , 
vary  one  geometric  parameter  to  provide  more  information  for  future 
designs  and  test  the  device. 

1.3  The,  1-ultiele  Inlet  Vortex  Valve 

A  multiple  inlet  vortex  valve  is  a  cylindrical  cavity  with  two 
opposing  sets  of  tangential  inlets  arranged  so  that  with  either  set 
flovjing  alone  maximum  vortex  strength  is  obtained.  Figure  1.0 
illustrates  this  geometry.  The  clockwise  inlets  are  suppled  from 
a  plenun  above  while  the  counter-clockwise  inlets  are  supplied 
from  a  plenum  below.  The  int ere  "•ting  feature  provided  by  the 
geometry  is  its  input  impedance.  Input  impedance  is  defined  as 
the  slope  of  the  curve  of  pressure  versus  mass  flow  at  the  counter- 
clockwise  ports  while  the  clockwise  ports  are  supplied  with  constant 
mass  flow  through  a  choked  orifice.  The  resulting  characteristic  is 
as  shown  in  Figure  1.1  Positive  flow  is  defined  as  flow  into  the 
vortex  cavity. 
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A  discussion  of  this  characteristic  wiH  lead  to  a  batter 
unders  banding  or  the  operation  of  the  valve .  Suppose  that  tl 
flow  set  of  inlets  cause  a  clockwise  vortex  while  flowing  alone 
therefore  the  other  inlet  set  would  provide  a  counter-clo  . 
vortex  while  flowing  alone.  For  discussion  purposes,  call  the 
counter-clockwise  inlets  number  1  and  the  clockwise  inlets  number 
2.  Referring  to  Figure  1#1,  region  1  corresponds  to  flow  out  of 
the  vortex  cavity  through  the  number  3.  inlets.  Point  "B"  corre- 
sponds to  the  pressure  caused  by  the  flow  through  the  number  2 
inlets  with  no  flow  through  the  number  1  inlets.  As  flow  is 
introduced  into  the  number  1  inlets ,  it  decreases  the  vortex 
strength  and  causes  the  pressure  to  decrease;  thus,  region  2  of 
the  characteristic  is  generated.  At  point  "A"  the  angular  momenta 
of  the  number  1  and  number  2  inlets  canrrel.  As  more  flow  is 
introduced  into  the  number  1  inlets,  a  counter- clockwise  vortex 
is  formed  and  the  pressure  increases",  thus,  region  3  of  the 
characteristic  is  generated. 

This  characteristic  provides  a  negative  resistance  region  and, 
therefore,  a  relaxation  oscillator  could  be  built.  A  flow  source, 
a  horizontal  line  on  the  coordinates  of  Figure  1.1,  would  have 
stable  intersections  in  all  regions  of  the  characteristic.  How- 
ever, a  pressure  source,  a  vertical  line  on  the  coordinates  of 
Figure  1.1,  would  have  stable  intersections  in  regions  1  and  3  and 
unstable  in  region  2.  Suppose  there  is  a  disturbance  in  the  f3 
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source  which  causes  a  f low  increase  momentarj  3y  while  the  pressure 
remains  fixed,  in  region  2  the  pressure  required  by  the  number  1 

inlets  decreases  and,  therefore,  the  source  provides  more  flow 
until  a  stable  intersection  is  reached  vertically  above  in  region 
3. 

A  relaxation  oscillator  can  be  constructed  if  the  number  2 
in3.ets  are  supplied  through  a  choked  valve  with  a  constant  bias 
flow  while  the  number  1  inlets  are  supplied  through  a  volume. 
The  flow  to  the  volume  is  set  by  a  choked  valve  at  approximately 
one-half  the  bias  flov;.  The  volume  acts  as  a  pressure  source. 
This  combination  provides  excessive  flow  for  region  1  and 
insufficient  flow  for  region  3»  The  result  is  a  Unit  cycle  in 
which  the  volume  discharges  in  region  3  until  the  vortex  valve 
switches  to  region  1  where  the  volume  charges  until  the  vortex 
valve  switches  back  to  region  J,   and  the  cycle  repeats.  Figure 
1.2  is  a  schematic  diagram  of  the  relaxation  oscillator. 

1A    Scosg. 

The  scope  of  this  thesis  is  to  design  and  construct  a  proto- 
type to  examine  one  parameter  of  the  internal  geometry  of  the 
valve,  to  predict  and  measure  the  pressure  sensitivity  of  the 
oscillator ,  and  to  predict  and  measure  the  temperature  s - . 
ity  of  the  oscillator. 


-7- 


Win^0.5YV.  . 
in  Diac 


VOLUME 


W 
^bias 


OUTLET 


FIGURE    1.2     SCHEMATIC    OF    RELAXATION   OSCILLATOR 


-8- 

CHAPTER  II 
APPARATUS 

2.1  Valve  Design 

Most  of  the  design  lias  been  done  empirically  clue  to  the 
absence  of  analytical  methods.  The  valve  has  been  designed  so 
that  the  ratio  of  inlet  area  to  outlet  area  could- be  varied. 
Previous  designs,  see  Reference  (1),  were  constructed  so  that  the 
opposing  tangenti?.!  inlet  sots  were  l2"ing  in  the  sarie  plane. 
Therefore,  the  present  design  separated  the  sets  of  tangential 
inlets  axially  by  0.25  inches.  Both  of  the  geometric  variations 
xrill  provide  useful  empirical  information  for  further  design  and 
possibly  lead  to  a  correlation  for  analytical  design. 

The  vortex  chancer  diameter  was  chesen  to  be  2.0  inches. 
The  vortex  cavity  has  two  outlets',  each  outlet  area  is  0.109  so. 
inches.  To  provide  versatility  in  the  inlet  configuration, 
twenty-four  inlet  slots  were  provided  every   fifteen  degrees 
around  the  vortex  chamber  perimeter.  In  order  to  provide  twenty- 
four  inlets  aroung  the  circumference  of  a  two  inch  diameter  circle 
and  have  part  of  the  circumference  remaining,  the  inlets  must  mate 
an  angle  of  greater  than  35  degrees  with  a  tangent  to  the  circum- 
ference of  the  circle.  The  0.125  X  0.125  inch  inlets  were  machined 
at  a  ^-0  degree  angle  to  a  tangent  at  the  vortex  chamber  perimeter « 
Figure  2.0  is  a  detailed  drawing  shewing  only  ^  inlets  for  clarity. 
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The  inlets  are  supplied  from  plenums  above  and  below  part  2. 
The  volume  was  designed  as  an  integral  part  of  the  valve. 

This  design  assured  equality  of  the  volume  pressure  and  the  pressure 
at  inlets  number  1.  If  the  volume  is  chosen  very  large  compared  to 
the  vortex  cavity  volume,  the  amplitude  of  oscillation  corresponds 
closely  to  the  knees  of  the  curve  shown  in  Figure  1.1.  The  reason 
for  this  is  that  the  charging  tine  of  the  volume  is  much  greater 
than  the  s^jitching  time  of  the  valve.  As  the  volume  approaches 
the  volume  of  the  vortex  cavity,  the  oscillation  stops*,  presumably 
due  to  the  increasing  influence  of  the  valve  dynamics.  The  present 
device  T.jas  designed  with  the  volume  nearly  k   times  larger  than  the 
volume  of  vortex  chamber  cavity  to  assure  oscillation.  The  volume  = 
5.85  cu.  in. 

Figures  2.1  and  2.2  show  a  photograph  of  both  sides  of  part  1 
in  Figure  2.0. 
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FIGURE    2.1   BIAS    SIDE  OF  VORTEX  VALVE 


FIGURE   2.2    VOLUME  SIDE  OF  VORTEX  VALVE 
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2.2.  Test  Apparat1' 

Figure  2.3  is  a  schematic  of  the  experimental  setup.  Figure 
2,h  is  a  photograph  of  the  apparatus.  Helicoidal  test  gages  were 
used  to  measure  pressure  and  Fischer-Porter  variable  area  meters 
were  used  to  measure  flow  rates.  The  air  heater  was  a  k*5\?A 
Wiegand  gas  heater.  Temperatures  were  measured  with  copper- 
constantan  thermocouples.  Pressure  oscillation  was  measured  using 
a  Eynisco  PT  25-25  pressure  transducer.  Frequency  was  measured  by 
stopping  a  lissajous  figure  which  was  generated  with  the  transducer 
output  on  the  vertical  axis  and  a  sine  function  generator  output  on 
the  horizontal  axis  of  an  oscilloscope.  The  frequency  was  counted 
using  a  Hewlett-Packard  frequency  counter. 
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CHAPTER  ] 
EXPERIMENTAL  DETEj  I  iTION  0?  STATIC  VALY  !  C    CTERISTICS 

3  •  1  Experimental  Procejiure 

Referring  to  Figure  2.3 .  the  static  characteristics  of  the 
vortex  valve  were  obtained  as  follows: 

a)  The  flow  metering  pressure  was  set  and  held 
constant  at  ^-0  psig. 

b)  Bias  flow,"C" ,  was  set  and  held  constant  at 
either  11.1  SGPM  or  5.25  SCFK. 

c)  Region  1  of  the  characteristic  curve  was  obtained 
by  opening  valve  "3"  and  i^ecording  flow  at 
flowmeter  "B"  while  valve  ''A:|  remained  closed. 

d)  Regions  2  and  3  of  the  characteristic  curve 
were  obtained  by  opening  valve  "A"  and  record- 
ing flow  at  flowmeter  "A"  while  valve  B 
remained  closed. 

e)  The  pressure  at  the  number  1  inlets  was 
measured  with  a  mercury  manometer. 

The  above  procedure  was  repeated  for  each  of  the  test  conditions 
shown  in  Table  3»1» 

Leakage  flow  was  1.6% of  a  total  input  of  28.25  3C   , 
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3.2  Data 

Figures  3<-l  to  3«S  show  the  static  valve  characteristics  for 
the  test  conditions  of  Table  3.1. 


TABIE  3.1 

TS 

3T 

coKDrnor ; 

Test  No, 

hA. 

1  a  -•  s 

Flow  (S0?I^ 

Air*  Tcr.rarature  (c ."'") 

T  1 

0.30 

5.25 

65 

T  2 

0.30 

5.25 

250 

T  3 

0A5 

5.25 

65 

T  k 

0,^5 

5.25 

250 

T  5 

0.60 

5.25 

65 

T  i 

0.60 

5.25 

250 

T  7 

0.91 

5.25 

65 

T  8 

0.91 

5.25 

250 

T  9 

0.30 

11.10 

65 

T10 

0.30 

11.10 

250 

Til 

QM 

11.10 

65 

T12 

o.ks 

11.10 

250 

T13 

0.60 

11.10 

65 

Tl^l- 

0.60 

11.10 

250 

T15 

0.91 

11.10 

65 

T16 

0.91 

11.10 

250 
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3»3  -Analysis  for  Static  Gl  es_ 

The  analysis  for  this  device  was  originated  by  Ir.  Adam  C. 
Bell,  Reference  (l).  To  provide  continuity,  to  this  report,  it  is 
repeated  here* 

See  Figure  3»9»  The  governing  equations  for  the  device  i 
as  follows! 

a)  Source  Equation 

ty*   =  constant  (l) 

b)  Volume  Equation  (continuity) 

W.      -  W,      =   />V     do  (2) 

IF    dt 

c )  Load 

V^     *    g  (p^  T-   ;         .try)  (3) 

The  volume  is  an  integral  part  of  the  valve ;  therefore , 
p  =  p.  Since  IT-,  cannot  be  determined  as  an  analytic  function  of 
pressure,  p,  the  equation  (2)  cannot  be  integrated  directly. 
Therefore,  a  graphical  technique  will  be  used  for  the  solution. 
The  technique  is  slopa-line  integration  -  see  Appendix  3. 

Rewriting  the  volume  equation  in  finite  difference  form! 

6P  =  (win  -  w^.  sr  w 

where        6T  =  _B_  ot  (5) 
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This  form  of  the  equation  states:  "the  change  in  pressure  over 
an  increment  of  time  is  equal  to  the  average  not  flow  into  t 
chamber  during  that  tine  multiplied  by  an  integrating  tJU  a  fac  -  ". 
The  time  constant  of  the  volume  is  given  by: 


t  =  pJL  (6) 

B 


Frequency  of  oscillation  is  inversely  proportion'.'.],  to  tl 
time  constant  of  the  volume.  Since 

B  -  kp  (?) 

and  £  =  L  (8) 

P    RT 

the  time  constant  of  the  volume  is  inversely  proportional  to 

temperature.  If  the  static  characteristics  were  independent  of 
temperature,  we  could  predict  that  frequency  would  be  directly 
proportional  to  temperature.  Since  the  static  characteristics 
do  vary  with  temperature,  the  temperature  sensitivity  prediction 
becomes  more  involved.  Furthermore,  the  frequency  is  insensitive 
to  pressure.  However,  this  is  only  true  if  mass  flows  are  held 
constant. 
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3»^  Ijynajaic  3el    r»  T   'otical 

The  predicted  temperature  sensitivity  was  obtained  by  assu-.- 
ing  a  linear  variation  between  65°   F  and  250°  Fe  The  frequency  at 
the  end  points  was  predicted  using  the  graphical  technique 
mentioned  in  Section  3»3  and  illustrated  in  Figure  3»10»  Frequsn- 
is  calculated  by  relating  the  ST  period  of  oscillation  to  real 
tine  through  equation  (5)»  Pressure  history  can  also  be  predict 
by  plotting  pressure  versus  • 

For  purposes  of  illustrating  the  method,  the  predicted  fre- 
quency of  oscillation  for  Figure  3*10  is  calculated.  Figure  3»1° 
shows  the  period  of  oscillation  to  be  25  ST  •  The  slope  of  the 
lines  is: 

ST  =  2.0  psig  (  9) 

SCFM 

Relating  this  period  of  oscillation  to  real  time 

St     =    JL     ST  (10) 

kRT 

where  V  =  3-385  x  10"3  cu.  ft. 

k  =  1.4 

T  n  525  °R 

R  =  1.72  x  103  ft. 2/  sec. 2  OR 
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6t  =  bS»5   m  sec.  therefore 

f  =  20,2  Hz 

The  results  of  the  integration  of  the. remaining  stotic  character- 
istic curves  are  shown  in  Figure  3»H« 
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CHAPTE3  IV 
EXPERIMENTAL  RESULTS 

4.1  Geometric  Effects  • 

The  geometry  of  the  valve  has  many  effects  on  its  static 
characteristics.  The  static  characteristic  changes  considerably 
as  the  inlet  to  outlet  area  ratio  increases.  The  characteristic 
also  changes  with  temperature.  As  the  characteristic  changes, 
so  does  the  frequency  and  the  output  waveform.  See  Figures  4.1 
and  4.2.  Furthermore,  the  temperature  sensitivity  of  the  device 
can  be  changed  by  changing  valve  geometry  as  shown  in  Figure 
3.10.  The  geometric  configuration  where  A./A0  =  0.60  predicts 
a  34.6  percent  change  in  frequency  for  a  35*2  percent  change  in 
absolute  temperature . 

The  position  of  the  inlet  slots  around  the  vortex  chamber 
were  varied  holding  A^/Aq  constant.  It  was  found  that  spacing 
the  inlets  symmetrically  or  assymmetrically  had  little  effect 
on  the  valve  characteristic,  provided  the  two  sets  of  tangential 
inlets  remained  axially  opposed. 

The  frequency  of  oscillation  should  be  inversely  proportion- 
al to  the  volume  for  a  given  vortex  cavity  geometry.  Scaling  down 
of  the  entire  device  should  increase  frequency  if  the  laws  of 
dynamic  similitude  apply.  However,  it  is  not  clear  how  scaling 
down  the  vortex  cavity  effects  the  static  characteristic  of  the 
device . 
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FIGURE    4.1     VOLUME    PRESSURE    HISTORY 
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^•2  Iressure  Sensitj 

The  device  is  very  nearly  insensitive  to  supply  pressure  if 
air  temperature  remains  constant.  The  actual  pressure 
of  the  device  is  shown  in  Figure  4.3*  The  sensitivity  was  Measured 
by  setting  flow  conditions  at  a  metering  pressure  of  40  psig. 
Supply  pressure  was  then  increased  to  80  psig,  then  docre-:.:;ed  to 
10  psig.  The  ratio  of  flow  into  the  volume  to  the  bias  flow 
remained  nearly  constant  tliroughout  the  range  tested. 

4.3  Temperature  Sensitivity 

The  ratio  of  predicted  temperature  sensitivity  to  actual 
temperature  sensitivity  was  3»2  to  1  for  a  given  geometry*  The 
main  problem  encountered  while  attempting  to  measure  temperature 
sensitivity  was'  a  controllable  constant  temperature  source.  The 
fact  that  such  apparatus  was  not  available  made  it  impossible  to 
obtain  steady  state  temperature  data.  A  further  shortcoming  of 
the  apparatus  was  that  the  device  was  constructed  of  brass  and, 
therefore,  provided  a  large  heat  sink. 

However,  temperature  sensitivity  data  was  obtained  while  the 
air  temperature  was  steadily  increasing.  Though  this  is  highly 
undesirable,  it  was  the  best  information  capable  of  being 
from  this  apparatus.  To  guard  against  misrepres     ;:  the  temp- 
era tore  sensitivity;  temperature  was  measured  Lnlet, 
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valve  outlet,  and  in  the  large  volume.  Table  k»1  cl 
the  large  temperature  gradients  present. 

Figure  k-J\-   shows  the  actual  temperature  sensitivity  of  the 
device,  where  temperature  of  the  air  in  the  volurse  is  plotted. 
Although  the  sensitivity  does  not  agree  with  that  predicted  in 
Section  3«^i  "the  author  feels  that  the  undesirable  thermal 
characteristics  of  the  valve  are  a  major  cause  of  deviation. 
The  predicted  gain  was  obtained  with  the  entire  device  at  either 
65°  F  or  250°  F.  However,  steady  state  temperatures  were  not 
obtained  during  the  temperature  sensitivity  tests.  Further 
pursuit  of  temperature  sensitivJ.ty  data  with  this  apparatus  v:as 
not  attempted. 
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TABUS  4.1 

air  temp 

Inlet  Temp  (°?) 

Outlet  ?er:v   (°F) 

Volume  Teino  (°F) 

65 

6A 

65 

80 

77 

70 

9? 

92 

80 

326 

106 

90 

126 

117 

100 

139 

127 

110 

152 

140 

120 

162 

150 

130 

171 

159 

140 

181 

170 

150 

191 

181 

160 

200 

190 

170 

210 

201 

180 

219 

211 

190 

226 

220 
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240 
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CHAj 
CONCLUSIONS  A?JD  SUGGESTIONS 

Geometric  parameters  have  a  significant  influence  on  the  va"1 
characteristics .  Herein  li.es  an  area  where  much  more  investigation 
should  be  conducted.  The  effect  of  each  of  the  geometric  parameters 
should  be  studied  to  lead  eventually  to  an  optimum  design.  The 
effect  of  scaling  the  device  is  vitally  important  to  its  eventual 
application  and,  therefore,  should  be  studied* 

This  work  has  shown  the  device  to  be  in  fact  temperature 
sensitive  and  pressure  insensitive.  Future  work  should  begin  by 
constructing  a  mold  from  the  existing  device  and  casting  a  r 
valve  using  a  non-heat-conducting  material,  possibly  plaster. 
Prior  to  attempting  to  measure  the  temperature  sensitivity,  a 
controllable  temperature  source  must  be  obtained.  If  the  results 
are  acceptable,  a  device  should  be  constructed  of  a  material  which 
follow  the  guidelines  established  in  Appendix  A.  If  a  desirable 
temperature  sensitive  device  is  designed,  it  must  be  tested  for 
speed  of  response  to  step  changes  in  temperature. 
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AP]     .  A 
HEAT  TRANS  ] 

To  be.  a  useful  temperature  sensor,  it  is  desirable  to  sense 
the  temperature  of  the  air  entering  the  device.  The  following  is 
an  order  of  magnitude  study  to  determine  a  material  which  results 
in  the  least  heat  transfer  to  the  device. 

For  this  study,  assume  a  solid  homogenous  body  whose  physical 
properties  are  constant;  it  is  initially  at  a  uniform  temperature 
Tjl»  and  is  suddenly  plunged  into  a  fluid  whose  temperature  is 
uniform  at  T~;  the  heat  transfer  coefficient  at  the  surface  of 
the  body  remains  constant. 

Assume  that  the  Biot  number  (k/nr0)  is  greater  than  6.  With 
this  assumption,  the  internal  temperature  gradients  are  negligib"1^, 
i.  e.  the  plate  is  thin. 

An  energy  balance  may  be  written  for  the  body: 

Heat  input  =  increase  of  intern?.!  energy 

dQ  =  Ah  (Tf~T)  dt  b  cyoV  dT  (A.l) 

With  the  assumptions  made,  integration  of  (A.l)  yields* 

t  =  cpV     In  Tf  -  Tx  (-.2) 

•   hA      Tf  -  T2 

Equation  (A. 2)  gives  the  time  required  to  heat  the  body  f . 
temperature  T-^  to  T2« 
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Solving  equation  (A.2)  for  the  time  to  rise  ( . 
the  way  between  an  initial  and  final  temperature  yields  the  tin 
constant: 

t  =  /ocV  (A.3) 

For  the  temperature  sensor,  we  want  the  material  tc  :    :>nd 
to  a  change  in  temperature  quickly* 

For  flowing  gases,  the  order  of  magnitude  for  h  is  2  -50. 
The  actual  value  i-rill  depend  on  the  gas.  The  time  constant  decreas- 
es when  the  device  is  scaled  down  since  V/A  decreases  with  scaling. 
The  best  material  to  use  to  construct  the  device  will  have  the 
smallest  /oc  product  with  the  constraint  that  the  Biot  number  is 
greater  than  6.  Table  A.l  shows  the  thermophysical  properties  ot 
1112°  F.  The  half  thickness,  r0,  of  the  material  for  this  type 
device  is  on  the  order  of  0,02  ft.  Assuming  h  to  be  50?  the 
thermal  conductivity,  k,  for  a  Biot  number  equal  to  6  is  on  the 
order  of  6.  The  po   product  for  the  various  alloys  listed  in  Table 
A.l  are  the  same  order  of  magnitude. 

The  design  of  an  operational  device  could  isolate  the  device 
from  its  environment  by  constructing  a  plenum  chamber  around  the 
entire  device  and  discharge  the  exhaust  gas  from  the  valve  into 
the  plenum.  This  design  should  greatly  reduce  t     ature 
gradients  in  the  device. 
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TABIE  A.l 
THEKMOFHISICAL  PROPERTIES  AT  11120  F 

/oc  (DTU/ft.3  ojr)  k  (BTU/hr.ft.°F) 

Copper,  pure             51.1  20^ 

Iron,  pure              53*3  23 

Molybdenum              38.2  6l 

Nickel                59.0  32 

Steel,  mild             55.0  19 

Stainless  steel           53  »5  13 

Tungsten                38.6  65 
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The  author  fe.^ls  a  device  constructed  o^  stainless  steel 
(k  -  13)  and  designed  to  be  isolated  from  its  environment  should 
minimize  heat  transfer  to  the  device  and  thereby  permit  the  devj 
to  be  used  as  a  sensor  of  input  air  temperature* 
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APPENDIX  B 
SLOPS-LINE  INTEGRATION 

This  appendix  describes  slope-lino  integration  in  gener< '. 
and.  explains  its  use  in  analyzing  the  static  valve  characteristics. 

Reference  (6)  is  the  source  of  the  choral  description  of  the 
technique . 

See  Figure  B.  1.  Suppose  we  are  given  tx;o  rectangles,  A, 
and  A£,  having  equal  bases  Ax.  The  areas  of  the  rectangles  are 
therefore  proportional  to  their  heights,  h,  and  h?. 

h-.  e  !n  (b.i) 

A2    h2 

Suppose  v;e  are  also  given  two  inverted  isosceles  triangles 
each  having  the  sans  vertex  angle.  The  bases,  a-,  and  a?,  of  the 
inverted  isosce3.es  triangles  are  also  proportional  to  their 
heights,  h-i  and  h£» 


*1  ~  \  (B.2) 

a2     h2 


Therefore ,  the  bases  of  the  inverted  isosceles  triangles  are 
proportional  to  the  areas  of  the  corresponding  rectangles 

l2' 


A,  and  A0. 


■':'; 


K, 


_JL 


A, 


K 


h-3r  H 


K— /^,  X H  f<— Z\  X 


FIGURE      B. 


A, 


A- 


K— ^  X  -^  ^— ^X 


FIGURE      B.2 
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t3L   =  -i  =  !l  (B.3) 

A2     h2     a2 


Suppose  A^  and  A2  were  superimposed  as  shown  in  Figure  3.2. 
They  are  in  effect  inscribed  in  and  circumscribed  about  a  trape- 
zoid Ay     The  area  of  A~  =  0.5  (h^  +  h2)Ax  which  is  equal  to  the 
average  area  of  the  two  rectangles  A-  and  A?;  A,  =  0.5  (A-  +  A?). 
The  area  Aq  can  be  found  graphically;  therefore,  by  adding  half 
of  the  "V"  used  to  determine  A-  to  half  of  the  "V!l  used  to  deter- 
mine A2»  The  result  is  a  "V"  with  a  linear  distance  ao  which  is 
a  measure  of  area  A«.  See  Figure  B.2. 

A  curve  can  be  approximated  by  a  series  of  trapezoids; 
therefore,  the  area  under  the  curve  can  bo  found  graphically.  If 
the  widths  of  the  trapezoids  are  all  the  same ,  a  series  of  V's 
all  having  equal  vertex  angles  will  determine  the  area  of  the 
trapezoids.  See  Figure  B#3«  Note  that  the  points  in  the  slope- 
line  diagram  record  the  growth  of  area  plotted  against  values  of 
y,  not  values  of  x. 

Suppose  we  assume  that  the  procedure  is  exact.  The  assump- 
tion is  valid  if  the  curve  is  a  polygon  rather  than  a  smooth 
curve.  The  procedure'  is  more  precise  for  a  smooth  curve  as  the 
width  of  the  trapezoids  approaches  zero,  i.  c.  steep  slopes. 
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Since  dA  =  ydx  (3.^) 

y  =  dk  .5) 

dx 


Therefore),  Figure  B*3  can  be  viewed  as ' 

A'  =  dA  =  f  (A)  (B.6) 

dx 

Such  a  graph  is  a  phase-plane  plot. 

Now,  let  us  suppose  we  were  given  a  graph  of  A'  as  a  function 
of  A  and  we  wish  to  obtain  A  as  a  function  of  x.  We  first  choose 
a  Ax  interval.  We  then  can  set  up  the  two  slope-lines  that 
constitute  the  V,  In  Figure  B.^  we  see  that  AA  =  A'  Ax  since  the 
are--'  of  a  rectangle  A'  Ax  can  be  represented  by  AA.  We  now  start 
at  any  point  on  the  curve  and  construct  the  slope-lines  from  the 
given  curve  to  the  axis  and  back  to  the  curve,  etc.  We  obtain  a 
series  of  points  on  the  curve  spaced  Ax  units  apart.  T.\Te  can  now 
calibrate  these  points  with  appropriate  values  of  x.  When  we  have 
finished,  we  can  plot  A  as  a  function  of  x.  We  have  successfully 
solved  the  differential  equation  (B.6).  We  can  start  anywhere 
along  the  curve  since  this  reflects  the  arbitrary  choice  of  the 
constant  of  integration. 

The  direct  application  to  the  analysis  for  the  static 
characteristic  is  discussed  in  the  following  paragrc  '  ■  The 
analogy  that  can  be  drawn  between  the  present  analysis  and  t] 
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previous  discussion  is  straightforward.  The  static  valve  c1 

teristics  together  with  the  flow  into  the  volume,  V,Tin,  r     snt 

a  phase-plane  plot.  Since  the  not  flow  rate  into  the  volume, 

Win  "  W3  =T^P  is  plotted  against  p.  Therefore,  the  analogous 
dt 

quantities  are  I 


A 

P 

X 

t 

dA 

d  : 

dp 

dt 

Sj.nce  we  are  given  a  graph  of  dp  as  a  function  of  p,  v;e  can  plot 

dt 

p  as  a  function  of  t  and,  therefore,  can  obtain  the  frequency  of 

oscillation. 

The  procedure  is  to  establish  slope-lines  of  slope  8T/2  as 
shown  in  Figure  3»10  and  construct  the  zig  zag  pattern  between 
W^n  and  Wn.  The  linear  distance  between  the  V's  represent  6p. 
We  obtain  a  scries  of  points  on  the  curve  spaced  6T  units  apart* 
We  can  now  plot  p  as  a  function  cf Tand  thus  determine  the 
period  of  oscillation.  The  period  can  be  related  to  real  t 
through  the  equation 

St  =  JL  ST 
kRT 


and  thus  the  frequency  of  oscillation  can  be  computed. 
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